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1 Abstract 

In this article we review existing literature on dynamic copulas and then propose an n-copula which 
varies in time and space. Our approach makes use of stochastic differential equations, and gives rise 
to a dynamic copula which is able to capture the dependence between multiple Markov diffusion 
processes. This model is suitable for pricing basket derivatives in finance and may also be applicable 
to other areas such as bioinformatics and environmental science. 

2 Introduction and motivation 

Mapping joint probability distribution functions to copula functions is straight forward when they are 
static, due to Sklar's Theorem. On the other hand, mapping time dependent probability functions 
to copula functions is more problematic. In this article we 

(a) review the techniques for creating time dependent copulas and 

(b) extend the method described in [4], [5], since it incorporates both time and space. These 

equations are the first of their kind in higher dimensions, since only 2-dimensional examples 
have previously been described. 

There are at least two areas in which the time dependent copulas of this chapter are applicable, 

• Credit derivatives. We would assume in this application that we have a portfolio of n firms, 
and Xi{t) is the value of a i-th firm's assets at time t. Each marginal distribution associated 
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with Xi{t) would represent the probabihty of the firm's vahie falHng below some threshold, 
given certain information at time zero. The time varying copula would represent the evolution 
of the joint distribution or state of the entire portfolio. 

• Genetic drift. For example, each Xi{t) may represent the frequency of a particular gene at 
time t. Each marginal distribution would represent the probability that the frequency of a 
particular gene had fallen below some threshold. The copula would relate to the evolution of 
a group of genes of interest. 

2.1 Notation and Definitions 

In order to understand some of the issues surrounding the mapping of copulae to distributions it is 
necessary to go back to some of the basic definitions and some notation in relation to the probability 
distributions of interest. 

The notation used for a univariate probability transition function here will be 

Pr{X,(i,) < I Xj{tj) ^ Xj} = F{t^,x^ I tj.Xj). (1) 

If tj = then it is quite common to suppress the zero and so the notation the distribution in this 
case would be F{ti,Xi \ x^). 

2.2 Method of Darsow et al 

Authors in [I] were the first to attempt to map a transition probability function to a copula. First 
let us recall the definition of a bivariate copula. 

Definition 1. 2-copula. A function C : [0, 1]^ — > [0, 1] is a copula if it satisfies the following 
properties; 

1. C(ui,0) = 0, C(0,U2) = 0, for aU ui,U2 £ [0,1] 

2. C(ui, 1) = ui, C(l, U2) = U2, for all wi, U2 G [0, 1] and 

3. For every Ua, Uf,, Va, Vf, e [0,1], such that Ua < Ub, Va < Vb, the volume of C, Vc'([ua,Wb] x 
[va, Vb]) > 0, that is 

C{ub, Vb) - C{ub, Va) ~ C{Ua, Vb) + C{Ua, Va) > 0. 
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Sklar's Theorem. Suppose H is a bivariate joint distribution with marginal distributions Fi and 
F2 then there exists a 2-copula C, such that for all xi,X2 G M 

Hixi,X2)^C{Fi{xi),F2{x2)). (2) 

If Fi and F2 are continuous distributions then C is unique, otherwise C is uniquely determined on 
RanFi x RanF2, see Nelsen [B]. 

Definition 5.1. Markov Property. A stochastic process Xi{t) and Xi ^ R, a < t < b is said to 
satisfy the Markov property if for any a < ti < t2 ■ . ■ , < tn < t, the equality 

Pr{X,it) <x,\X,{h),X,it2),...,X,it^)} = Pr{X,(i) < x, \ X,{t„)} 

holds for any 2;^ G R. A stochastic process is called a Markov Process if it satisfies the Markov 
property described in Definition 5.1. 

The following notation will be used for an unconditional probability function at time ti > 0, 

MMti) < ^ FuiXr) (3) 

for a stochastic process Xi{ti) and Xi G M. Let 

OXi 

then the corresponding density function / in this case is such that 

A univariate transition probability function F (Markov process) can be mapped to a bivariate copula 
C by setting 

F{U,x^ I t,,Xj) = Vu,C{Ft^ix^),Ft^{xj))) , (4) 

where is the partial derivative with respect to the second argument of C. This mapping enables 
us to build in time, see [T]. The first marginal distribution is associated with time ti and the second 
with time tj. We take the partial derivative of the copula, since the probability to which it is 
mapped is conditional. This method is particularly useful for building Markov chains. One of the 
most important innovations which enabled the authors in IJ to link copulas to Markov processes 
was to introduce the idea of a copula product; 
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Definition 5.2. Copula product. Let Ca and Cb be bivariate copulas, then the product of Ca and 
Cb is the function Ca * Cb : [0, 1]^ — * [0, 1], such that 

{Ca*Cb)(x,y) = f V,Ca{x,z)V,Cb{z,y)dz. (5) 

JQ 

This product is essentially the copula equivalent of the Chapman-Kolmogorov equation, as stated 
in Theorem 3.2 of J]. We restate that theorem here (with modified notation). 

Theorem 3.2. Let Xi{t), t E T he a real stochastic process, and and for each s,t G T let Cst denote 
the copula of the random variables Xi{s) and Xi{t). The following are equivalent: 

1. The transition probabilities F{t, A \ s,Xs) — Pr{Xi{t) G A | Xi{s) — Xg} oi the process satisfy 
the Chapman-Kolmogorov equations 

Fit, A\s,x,)^ [ Fit, A I u, OFiu, \ s, x,) (6) 

for all Borel sets A, for all s < < G T, for all w G {s,t) CiT and for almost all Xg G M. 

2. For all s,u,t G T satisfying s < u < t, 

Cst — Csu * Cut- (7) 

The work in [T] has advanced both the theory of copulas and techniques for building Markov pro- 
cesses. This method is also used in [9 to formulate a Markov chain model of the dependence in 
credit risk. The discrete stochastic variable Xiit) is interpreted as the rating grade of a firm at a 
particular point in time. A variety of copulas were fitted to the data and gave mixed results. There- 
fore, no copula was the best for all data sets. This type of mapping of the transition distribution 
to the copula is very simple, however, one consequence is that an n-dimensional transition function 
requires a 2n-dimensional copula. In other words, as the dimension of the copula increases, the 
calculation of the transition function becomes more and more computationally cumbersome. 

The method in [T] has also been extended in so that an n-dimensional Markov process can be 
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represented by a combination of bivariatc copulas and margins. Hence, 
Pr{X.(ii) <2:i,...,X,(t„) <x„} 

71 

= Y[Pr{X,{U) < X, I X,{ti) = xi, . . . ,X,{U^i) = x,^i}Pr{X,{ti) < xi} 

1=2 

TL 

= [|Pr{X,(t,) < X, I X,{U.i) = x,^i}Pr{X,{ti) < xi} 

i=1 

2.3 Conditional Copula of Patton 

Another approach to building time into a copula was formulated in f7]. In order to explain this 
approach, we need to recall more definitions and set up notation. Firstly, let be a filtration, then 

Vt{X,<x,\T} = F,{x,\T) (9) 

The multivariate analogue of equation ([S]) is 

Pr{X <x\T}^ H{x I T), (10) 

for X = (xi, X, . . . , XnY" such that the volume of iJ, Vh{R) > 0, for all rectangles R £ R" with their 
vertices in the domain of H , see [8], 



H{+oo, Xi, +00, . . . , +00 I T) = Fi{xi \ T), and 
H{—oo,Xi,...,Xn\T)—0 for all a;i, . . . , a:„ £ M. 

Here Fi is the z-th univariate marginal distribution of H . See [7] for a bivariatc version of H . As 
expected, the density of the conditional H is 

h{x I T) = V,,,...,,„ff(x I T). (11) 

In equation ([9]), the distribution is atypical since it may be conditional on a vector of variables, not 
just one, as opposed to a typical univariate transition distribution. 
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The author in f7^ mapped the conditional distribution H{x \ J^), defined above, to a copula of the 
same order. That is, for all G M and i = 1,2, ... ,n, 

H{xi,...,x„\T)^CiFiixi \T),F2{X2 \ T), . . . , F„{x„ \ T) (12) 

J- is a sub-algebra or in other words a conditioning set. Such conditioning is necessary for C to 
satisfy all the conditions of a conventional copula. The relationship between the conditional density 
h and copula density c is 

n 

h{xi,X2,...,Xn I J^) =c{ui,U2,...,Un \ J^)Y\_fii^i I > (13) 

1=1 

where Ui = Fi{xi \ J-), i = 1,2, ... ,n and fi, i = 1,2, ... ,n are univariate conditional densities. 

In terms of time varying distributions, we can think of the conditioning set as the history of all the 
variables in the distribution. In the case of the Markov processes, it is only the last time point which 
is of importance. The implication of this type of conditioning is that the marginal distributions 
in the copula can no longer be typical transition probabilities, but are atypical conditional proba- 
bilities. Hence, if each Xi represented the value of an asset at time t, the associated distribution 
Fi would represent the distribution of Xi, given that we knew the value of all the assets in the 
model, Xi, X2, . . . , Xn, at some previous time, for example t — 1. In other words, we can rewrite 
the time- varying version of the distribution and copula above as 

Htix],xl...,x'l \Tt-i) 

= Ct{Fl{x\ I Tt-i),F^{xl I Tt-i), . . .,F^{x^ \ Tt-^) \ Tt-i), (14) 

where 

J-t-l — (T{Xf._i, Xf._i, . . . , Xf._i, Xf._2, Xf_2, . . . ,Xf._2, . . . ,Xi,Xi, . . . ,Xi). 

In [7] , the marginal distributions are characterized by Autoregressive ( AR) and generalized autore- 
gressive conditional heteroskedasticity (GARCH) processes. Ultimately, they are handled in the 
same way as other time series processes. 

2.4 Pseudo-copulas of Fermanian and Wegkamp 

As we have seen above, Markov processes are only defined with respect to their own history, not 
the history of other processes. Therefore, the method in [7] is good for some applications but not 
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practical for others. If we want marginal distributions of processes, conditional on their own history, 
for example Markov processes, and want to use a mapping similar to that shown in [7], then it 
is possible via a conditional pseudo-copula. Authors in 3J introduced the notion of conditional 
pseudo-copula in order to cover a wider range of applications than the conditional copula in [7 . The 
definition of a pseudo-copula is 

Definition 5.3. Pseudo-copula. A function C : [0, 1]" [0, 1] is called an n-dimensional pseudo- 
copula if 

1. for every u £ [0, 1]", C(u) — when at least one coordinate of u is zero, 

2. C(l,l,...,l) = 1, and 

3. for every u, v G [0, 1]" such that u < v, the volume of C, Vc > 0. 

The pseudo-copula satisfies most of the conditions of a conventional copula except for C(l, 1, Uk, 1, . . . , 1) = 
Uk, so the marginal distributions of a pseudo-copula may not be uniform. The definition of a con- 
ditional pseudo-copula is 

Definition 5.4. Conditional pseudo-copula. Given a joint distribution H associated with Xi , X2, . . . , Xn, 

an n-dimensional conditional pseudo-copula with respect to sub-algebras 

!F = {^1,^2, . . . , J-n) and ^ is a random function C(- 1^,0) : [0, 1]" [0, 1] such that 

H{xuX2,...,Xn) = C{Fi{xi I Ti),F2{x2 I J^2),---,F.n{Xn \ Tn) \ T , G) (15) 

almost everywhere, for every {xi,X2, . . . , XnY" & K", see [2]. 
2.5 Galichon model 

More recently, a dynamic bivariate copula was used to correlate Markov diffusion processes, see [1], [S]. 
Unlike the previous models of time dependent copulas, this model addresses the issue of spatial as 
well as time dependence. The model uses a partial differential approach to obtain a representation 
of the time dependent copula. An outline of the main result follows. Consider two Markov diffusion 
processes Xi{t) and X2{t), t G [0,T], which represent two risky financial assets, for example options 
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with a maturity date T . The diffusions are such that 

dXi{t) = pLi{X{t))dt + di{X{t))dBi{t) 

dX2{t) = fl2{X{t))dt + d2{X{t))dB2{t) 

dBi{t)dB2{t) = pi2{Xi{t),X2{t))dt, (16) 

where X{t) — {Xi{t)^X2{t))'^ , /ii, (Ti, for i = 1,2, are the drift and diffusion coefficients, respec- 
tively. The Brownian motion terms are correlated with coefficient pi2 G [—1,1]. One would like 
an expression for the evolution of a copula between the distributions Fi, F2 of Xi{t) and X2{t), 
conditional on information at time t = Q, J-t^- Firstly a joint bivariate distribution H is mapped to 
a copula C, by 

H{t,xi,X2 I Tt^) = C(t,Fi(t,xi I Tt,),F2{t,X2 \ Tt„) \ Tt„) (17) 
then the Kolmogorov forward equation is used to obtain an expression for VfC. Letting 

ui = Fi{t,xi I Tto) and U2 = F2{t,X2 \ J^u,), Mi,W2 e [0,1], 

X = {xi,X2)'^ and shortening the notation for the copula to C{t,ui,U2), then the time dependent 
copula in 4] is 

VtC(i,Ul,U2) = ^al{x)fUt,Xi I J^to)^l^C{t,Ui,U2) + \dl{x)fl{t,X2 \ J^to)^l^C{t,Ul,U2) 

- VuiC{t,ui,U2)BiFi{t,xi \ Tt„) 



V„i,„2C(t,Mi,U2)/2(i, ^2 I J^to)BlFi{t, Zi I Tta)dz2 

( — 00,1:2] 

Vu^C{t,Ui,U2)B2F2{t,X2 I Tto) 



V„i,„2C(t,-Ui,U2)/l(i, Zl I ^to)B2F2it,Z2 I Tta)dzi 

' ( — 00,2:1] 

+ o-i(a;)CT2(x)pi2(xi,a;2)/i(i, xi | Tt„)f2it,X2 \ ^to)V„i,„2C(t, ui, W2), 



where Bi and B2 are the following operators, given any function g G C^(]R), 

Big = |v,,Qa2(a;)^ -;ii(a:)|v,,g+ Qa2(a;)^V2^g 
B2g = |v,2(^'?2(a^)) -M2(a:)| V,,g+ Qa2(a;)jv2^(7 



(18) 
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and 




For the greatest flexibility we would choose 



inf{xi : Fi{t,Xi \ Tto) > Ui} 



Fr^{t,u,\TuX e[0,l 



That is, ^ is the pseudo-inverse. If Xi{t) and X2{t) are individually Markov, that is, di and /i^ 
depend only on Xi, for i = 1, 2, then the formula for the time dependent copula simplifies to 

VtC(t,Ui,U2) = ^l[x)a2{x)pi2{xi,X2)fl{t,Xi \ Tto)f2{t,X2 \ Tta)Vui,U2C{t,Ui,U2) 



The main aim of this chapter is to extend some of these current models of time dependent copulas. 
We derive an n-dimensional version of the model in [T . A reformulation is also given, in which linear 
combinations of independent Brownian motion terms are used. 

3 n-dimensional Galichon Model for CDOs 

Suppose we have an n x n system of stochastic differential equations, such that X{t) e M" and B{t) 
is an n-dimensional Brownian motion. The vector X{t) could represent a portfolio of risky assets, 
as in a Collateralized Debt Obligation. We want to find a partial differential equation with respect 
to a time dependent n-copula, which gives us information on the riskiness of the package of assets. 
As in the 2-dimensional model, i is a scalar such that t G (0, T]. Let be a cr-algebra generated 
by {B{s);s < t} and assume X{t) is measurable with respect to J-t- In this case the diffusions are 
such that 



+ ^al{x)fUt,xi \ TtoM^C{t,ui,U2) + -al{x)fl{t,X2 \ TtoM^C{t,ui,U2). 



(19) 



dX{t) 



^{X{t))dt + AdB{t) 



(20) 



dB,{t)dBj{t) 



p^,{X,{t),X,{t))dt, 



(21) 
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where 



dX{t) 



f dXi{t)\ 
dX2{t) 



\dXn{t)J 



dB{t) 



a{X{t)) 



( dBi{t)\ 
dB2{t) 



\dBrXt)) 

l^i{X{t))\ 
a2{X{t)) 



VniX{t))J \^niXit))/ 

Note that in this case /i and a are n-vector functions which represent the drift and diffusion coeffi- 
cients of the process, respectively. Let A be 

faiiXit)) \ 

^2{X{t)) ... 



A = diag{a{X{t))) 



V 



The correlation coefficients py G 1] ^^^d let p be 

/ 1 pi2{Xi{t),X2it)) 

P2l{X2it),Xi{t)) I 

\pnl{X^it),Xiit)) 



. an{x{t)).J 
Pi„(Xi(i),x„(<))\ 

P2n{X2it),X„{t)) 

1 J 



Three conditions are required for the existence and uniqueness of a solution to equation (|20p : 

1. Coefficients p{x) and (j{x) must be defined for x e M" and measurable with respect to x. 

2. For x,y ^ R", there exists a constant K such that 

\\Kx)-Ky)\\ < K\\x-yl 
\\d{x)-a{y)\\ < if II a; -j; II, 
\\ f,{x) r + \\ a{x) r < K^il+Wxf) 
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and 

3. X{0) does not depend on B{t) and E[X(0)2] < oo. 

Theorem 5.1. The time dependent n-copula WiC{t,u) between a vector of distributions = 
Fi(t, Xi I .To), z = 1, . . . , n, associated with the Markov diffusions X{t) — \Xi{t)^ . . . , Xn{t)]'^ , condi- 
tional on information at time t = 0, Tt^ = xq is 

1 " /■ 

VtC{t, «) = 2 E / ^ii^yffi*' I a;o)V,„.,,„.,,„ V2^C(t, u)dz 

i=l •^(-oo,a:] 

+ ^ (- Vu,C{t,u)BiFi{t,Xi \xo)+ f V,,^.,i,,.,,Yu,C{t,u)BiFi{t,Zi \ xo)d 

i— 1 V J(-oo,x] 



1 " /■ 

ij=l^(-oo,i] 



(22) 



where Xi = ^{t,Ui \ Tto)i i = 1, . . . ,n. The intervals for the integration are 

(— oo, x] = (—00, xi] X . . . X (—DO, Xi-i] X (—00, Xi+i] X ... X (—00, Xn] and 
(— oo, a;] = (— oo, a;i] x . . . x (— oo, Xi-i] x (— oo, Xi+i] x . . . x (— oo, xj-i] 

X (-00, Xj+i] ... X (-00, Xn]. 

Also note that 

dz = dzidz2 ■ ■ ■ dzi-idzi+i . . . dzn-idzn and 

dz = dzidz2 ■ ■ -dzi-idzi^i . . .dzj-idzj^i . . .dZn-idZn- 

Thus, the i-th term is excluded in the first two integrals on the right hand side of Theorem 5.1. 
Similarly, in the last integral the i-th and j-th terms are excluded. Furthermore, for any smooth 
function g 

yZl,.,Zi,.,Zr,9 = 

and Vzi,.,zi,zj,.,zn9 = 



dz\ . . . dzi-idzi+i . . . dzn 



dzi . . . dzi-idzi+i . . . Zj^idzj+i ■ ■ ■ dzn 
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The operators Bl, i — 1, . . . ,n are the same as in the two dhxiensional model. If the diffusions 
are individuaUy Markov, that is, each ak and /ifc depends only on then the expression for VtC 
simplifies to 



1 " 



-Tr{ [n^{t, u) - diag{V2^C(t, u),\7l^C{t, i 



yl^Cit,u)}]DApA^D^}, 



where 



D = 



h 

h 










Inj 



Proof. In this case 1-dimensional Ito formula for each component of X{t) is 

dgiXS)) = {V..^g{Xm^l,{X{t)) + \v%g{X,{t))^UX{t)))dt 
+ V^^g{X,{t))a,{X{t))dB,(t). 

Define the vector Vx of partial derivatives with respect to components of x, as 

lv,,g{X{t))\ 



and the Hessian matrix of g{X{t)) 

n^ixit))^((v.,.^g{x{t)) 



l<'i,j<n 



(23) 



In this case, assume g £ C^(R"), then the n-dimensional Ito formula for g{X{t)) is 
dg{X{t)) 

= {(V,5(^W),/^(^W)) + ^Tr(H^(X(0)ipi^)}dt + V,5(^Wr^cfS(i), 

where (a, b) = a^b for any vectors a and b. Let the operators A on distributions (Kolmogorov back- 
ward equations), analogous to those in [4], [5], be called Al and for the 1- an n-dimensional case, 



12 



respectively. With respect to typical distributions Fi{t, Xi | r, and H{t, x \ r, the operators are 

I = n,{x)\/^Mt,Xi I T,^i) + ^a^VlF,{t,x, \ r,^,) (24) 

and 

A'^H{t,x I r,0 = {'7iH{t,x I T,0,^{x)) + ^Tr (nf{t,x \ T,^)ApA^) . (25) 

The operators Al, i = 1, . . . ,n and A^ are not used in the rest of the formulation, but are mentioned 
briefly, in view of the fact the Kolmogorov forward equations, which are required, are the associated 
adjoint operators of these. Assuming the density functions of H and F are h and /, respectively, 
then the adjoint operators A^* , i = 1, . . . , n and A^* have the form 

^*/i(i, Xi I T, ^i) = - Vx, [l^i{x)fi{t, Xi I T, ^i)] + [^a^fiit, Xi I T, ^i)] (26) 

and 

Arh{t,x\T,^) 

n 1 
i=l ^ij=l 

(27) 

The marginal density functions fi and joint density h, are such that 
fi{t,Xi I J^to) = fi{t,Xi I a;o), and H{t,x \ Tto) = H{t,x \ xq), where 

xq = {xi = Xi{Qi), X2 = -^^2(0), . . . ,Xn = Xn{0)), scc Appendix 5. A. In other words, the assumption 
made here is that all the distributions are conditional on the entire vector of realizations of x at 
time zero. As in the 2-dimensional case, it is possible to express the operator Af* in terms of the 
operators Al* associated with the univariate distributions; 

n 1 

A^*g = J2A*9 + 2 51 ^^-^i [pij{xi,Xj)ai{x)aj{x)g] . (28) 

i=l 

Given that 

ytfi{t,Xi I xo) = Ai*fi{t,Xi I xo), (29) 
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we can integrate the left hand side of with respect to Xi, call it B}, and we obtain 

BlF,{t,Xi \ Xq) = / \/tMt,Zi\ Xo)dzi 

J ( — 00,2;^] 



{ — 00,Xi\ 

VtF^{t,X^\xo). (30) 



Integrating the right hand side of with respect to Xi gives us 

/ A\* fi{t,Zi I xo)dzi = ~'^ii{x)fi{t,Xi I xo) + Va;i[^CTi(x)/i(t, I a;o)] 

J( — OO.Xi] ^ 



= [^xAIs-Ux)} - ^i^{x)]\/xMt,X^ I xo) + ^fff (x) V^^ (t, X, | Xq), 



Similarly, integrating over A^* will give us the analogous operator B]^ for the multivariate distribution 
H . Now, since 

2: I xo)dz 

(—00,2:] 

/ ^r^l^'-z I 2;o)rf2 + i ^ / Vz,,2^ [pjj(zi,Zj)d-i(z)CTj(z)/i(t, z I xo)](iz, 



■',3 



(32) 



where (— oo, x] = (—00, xi] x . . . x (— oo,x„], it is possible to get an expression for in terms of 
B\. That is, let Bl^H{t,x \ xq) = \7tH{t,x \ xq) and given that h{t,x \ xq) = V xi,.,x„H{t,x \ xo), 
we have 



"H{t,x\xQ) ^ 9 X! / '^z,,zj[pi3{zi,Zj)at{z)aj{z)\/zi....,z„H{t,z \ xo)]dz 

+ Al*V,,,...,,„H{t,z\xo)dz. (33) 

J{-cx>,x] 

The right hand side of equation (j33p can be expressed in terms in terms of the univariate operators 
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n „ 

+ By,,,.,,,,.^,^H{t,z\xo)dz. (34) 

^^-^ J{-oo,x] 

Let 

H{t,x I xo) = C{t,Fi(t,xi I xo),F2{t,X2 \ xq), . . . ,F„(t, a;„ | xo) \ xq) (35) 

where C is an n-copula defined on [0,T] x [0, 1]". At this point wc shorten the notation so that 
C{t,F{t,x I a;o)) is the same copula as above. We now seek an expression for B'^C{t,F{t,x \ xq)) 
by substituting for H with C in equation p4p . Letting Fi{t,Xi \ xq) — Ui, i — l,2,...,n, and 
u ~ (mi, . . . , UnY' , then from the first term in equation p4p we obtain (see overpage) 
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n „ 

SiV,„.,5„.,,„ff(i,z|xo)d^ 

i=l ''i-oo,x] 

n ^ 

i=l -'i-oo,x] 

+ I ^^lyz,,.,z,,.,z^C{t,F{t,z\xo))dz\ 

= Ef/ {Vz,4r^-Mi(2)}/i(i,^i ko)V„,V,„.,,„.,,„C(t,u)dz 

i=l VV(-<x,,x] ^ 
J(-oo,a;] ^ 

+ / 4r^V,Ji(t,0i |a;o)V„,V,„.,,„.,,„C(i,w)d^') 

J(-oo,S] ^ / 
^ Vi(-cx,,x] ^ 
J(-oo,ci;] ^ 

^(-oo,5] ^ / 
n ^ 

= E/ V.„.,5„.,.„V„,C(i,w)SjFi(t,^i |a;o)d^ 

8=1 -^(-00,2:] 

n „ 
i=l •^(-oo.x] 

Since z is a dummy variable and the multiple integrals exclude that over (—00, a;,], we can write 

n ^ 

8=1 -^(-OOiXl 

= E/ Vzi,.A,.,z„V„,C(i,u)BJFi(t,^i I a;o)rf^ 

j^;^ J(-oo,2:] 
1 " /■ 

+ oE / |a;o)V,„.,,„.,,„V^,C(t,«)rf^. (36) 

i=l -^(-oo.a:] 
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so 



From the second term in equation p4p we have 
1 " /■ 

^ ij=l-^(-oo,x] 
1 " /■ 

1 " /■ 

9 51 Ptj{Xi,X.j)d-^{z)aj{z)fi{t,Xi \xo)fj{t,Xj \xo)V z^^,^2i,z,,.,zyui,up{t,u)dz 



j^;^ J(-oo,2;] 
1 " /■ 
^ .-1 ^ f — oo.kI 



1 " /■ 

+ — ^ ^ Pij{Xi,X j)ai{z)aj{z)fi(t,Xi \xo)fj{t,Xj \xQ)Vzi,.,zi,z,,.,zyu,,up(t,u)dz. 



Now, we also have 



VtH{t,x\xo) = B^H{t,x\xo) 

n 

= ytC{t,F{t,x I xo)) +^ V„.C(t,F(i,a; | xo))Vti^,(t, | xq) 

n 

= \/tC{t,u) +Y,^u,C{t,u)BlF,{t,x., \ xo). 

i=l 
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Matching equation ([55)1 and ([5^ and rearranging, we obtain 
VtC{t,u) 

= J2(- ^u.Cit,u)BlF,{t,x, \xo)+ I V,,,.,i^,.,,^^u,Cit,u)BlF,{t,z, \ xo)dz 

i=l ^ J{—x>,x] 

1 " /■ 

1 " /■ 

If the equations are individuaUy Markov, so that each (7^ and /i^ depends only on Xk^ then 

^ {- Vu,C{t,u)BlF,{t,x, \xo)+ f V,,,.,,,,.,,„Vu,C{t,u)BlF,{t,z, \ xo)dz) = 0, 

so the expression for VjC simphfies to 
1 " 

i=l 

+ ^Tr{ u) - diag{V2^C(i, u),Vl^C{t, u),..., Vl^C{t, u)}]DApA^ D^} . 

□ 

4 Conclusion 

We have described a dynamic n-copula which varies in time and space. This copula is the first of its 
kind in greater than 2 dimensions. The dynamic 2-copula was previously described in [5]. In that 
case, the copula could be applied the pricing of pairs of options and other credit derivatives. In the 
n-dimensional case, it is possible to use the dynamic copula for the pricing of any basket derivatives 
or a number of commodities. Future work in this area may involve numerical experiments, sensitivity 
testing and simulations in order to determine how robust the model is. Other possible applications 
include that of the health industry and environmental science. 

I would like to acknowledge Alfred Galichon for his valuable discussions in relation to this work. 
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